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Abstract 

Using a sample of 3.3 x 10 6 i?-meson decays collected with the CLEO detector 
at the Cornell Electron Storage Ring, we have studied B~ — > D°£~v and 
B° — ► D + l~v decays, where l~ can be either e~ or \i~ . We distinguish 
B — ► ZWf from other I? semileptonic decays by examining the net momentum 
and energy of the particles recoiling against D — £ pairs. We find V(B — ► 
= (14. 1± 1.04: 1.2) ns _1 and derive branching fractions for B~ — > D®i~v 
and 5° -► D + l~v of (2.32±0.17±0.20)% and (2.20±0.16±0.19)% respectively, 
where the uncertainties are statistical and systematic. We also investigate the 
B — * Dlv form factor and the implication of the result for |V^,|. 
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The semileptonic decays of the P-meson play important roles in heavy quark physics. 
They provide our best information on the CKM matrix elements V c b and V u b and reveal 
the dynamics of heavy quark decay in their form factors. Heavy Quark Effective Theory 
(HQET) @ provides a framework for calculating the form factors for b — *> civ decays and 
suggests a reliable method for extracting \V c b\ by predicting that QCD effects are small at 
the kinematic point where the final-state meson is at rest relative to the initial meson (zero 
recoil). Most studies of the form factors and \ V c b\ [||] have used the B — > D*iv decay because 
its differential branching fraction near the zero-recoil region is large and the QCD effects 
have been calculated to the highest accuracy ||]. There have been two recent studies of the 
mode B — > Div || ||. Here we present a new, more precise study of B — > Div using a 
different analysis method. 

We select B charm semileptonic decays including B — > Div, B — > D*iv, B — > D**iv, 
and B — > D^ixiv decays by identifying events with a D (D° or D + and their charge 
conjugates) and a lepton. We then separate B — > Div from the other semileptonic modes 
using the net energy and momentum of the particle or particles recoiling against the D — I 
pair. Information on the partial width and differential decay rate are obtained from the 
B~ — > D°i~v and B° — > D + i~v yields, which we extract in bins of the HQET variable 
w = (Mg + Mp — 5 2 )/(2MbMd), where q 2 is the squared invariant mass of the virtual W. 

The data used in this analysis were accumulated in the CLEO detector J7| at the Cornell 
Electron Storage Ring (CESR). The data consist of 3.3 million BB events collected at the 
T(4f>) resonance. We also use 1.6 fb _1 of data collected 60 MeV below the T(4S') resonance 
to study the background from the e + e~ — > qq continuum. The CLEO detector measures 
the trajectories of charged particles in a drift chamber system inside a 1.5T superconduct- 
ing solenoid. The main drift chamber provides the specific ionization (dE/dx) of charged 
particles and their time-of-flight (TOF) is provided by scintillation counters surrounding the 
drift chamber. A Csl electromagnetic calorimeter is used in electron identification. Muons 
register hits in counters embedded in steel surrounding the magnet. 

In this analysis we select events having at least 5 charged tracks and, to suppress non- 
BB events, the ratio of Fox- Wolfram moments || H2/H0 < 0.45. We reconstruct D° and 
D + candidates in the decay modes D° — > K~tc + and D + — > K~h + 'k + respectively. We 
distinguish K's from 7r's based on a \ 2 probability (P# and P v ) that combines dE/dx and 
TOF. Each daughter track must satisfy Pj > 0.01 and Pj/(P^ + P n ) > 0.3, where % is n or 
K as appropriate. The candidate mass must be within 1.835 GeV < m KlT < 1.893 GeV for 
D° and 1.846 GeV < rriKmr < 1-890 GeV for D + decays. To suppress P-mesons produced 
in e + e~ cc events, we require |p_o| < 2.5 GeV/c. 

Electrons are identified using dE/dx, the shape of the shower in the Csl calorimeter, and 
E/p, the ratio of the candidate's energy deposit in the Csl to its momentum. Electrons are 
required to have momenta between 0.8 GeV/c and 2.4 GeV/c. Muons, selected within the 
same momentum window, must penetrate at least 5 interaction lengths of material. This 
requirement places an implicit lower bound on the muon momentum of about 1.4 GeV/c. 
We require that the lepton have the same charge as the K from the decay of the P-meson. 
Since the decaying B- meson is nearly at rest, the P-meson and the lepton are in opposite 
hemispheres for more than 90% of B — » Div decays; we demand that this be so. 

For events satisfying these criteria, we compute cos 9b-d£, the cosine of the angle between 
the D£ momentum p^e = Pd + Pe and the B momentum p#, assuming that the decay is 
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TABLE I. Yields and Backgrounds. Uncertainties are statistical only. 





D°£ 


D + £ 


Total Yield 


12595 ± 112 


18087 ± 134 


Rflnnnm f\Tr( it 1 ( ,nm nmntinnQ 






Uncorrelated 


948 ± 63 


761 ± 75 


Continuum 


452 ± 84 


432 ± 104 


Correlated 


119 ± 16 


119 ± 23 


Fake Lepton 


71 ± 19 


26 ± 25 


Background-subtracted Yield 


5922 ± 163 


3247 ± 201 



B — > Div and that the only missing particle is the massless neutrino; that is 

a 2E B E m -M 2 B -M 2 m 
cos9 B -Dt = 7Ti — rr-^i — ■ 1 

This quantity lies between —1.0 and 1.0 for B — > Div decays. When final-state particles 
are missing in addition to the u, as is the case for the other B semileptonic decay modes, 
cos 9b-d£ decreases. We use the distribution of cos9b-di to separate B — > Div events from 
the other B semileptonic decay modes after subtracting backgrounds. 

The backgrounds come from several sources: random Ktt^tt) combinations, D-mesons 
matched with a lepton from the other B decay (uncorrelated), D-mesons combined with a 
lepton that is a granddaughter of the same B (correlated), hadrons misidentified as leptons 
(fake lepton), and e + e~ — > qq events. 

Events in the mass regions above and below the D peak (sideband) are utilized to esti- 
mate the random combination contribution. A study using Monte Carlo-simulated B decays 
shows that this method subtracts the right amount of background within small uncertainties 
and that the cos 9 b-d£ distribution of the events in the sidebands reproduces that of the 
background events in the signal region. 

In uncorrelated background events, the angular distribution between the D and I is nearly 
uniform because they arise from different B- mesons, both nearly at rest. We take advantage 
of this uniformity: for each event in which the D and I are in the same hemisphere we reverse 
the lepton's direction and compute cos 9b-d£, thereby constructing this distribution for the 
opposite-hemisphere background events. 

The e + e~ — > qq continuum background is measured using events collected off resonance. 
The correlated background, which is small, arises from modes such as B — > D S D followed 
by D s — > Xlv and B — > DXtu followed by r — > Ivv. We estimate these contributions using 
a Monte Carlo simulation. Fake lepton background is estimated by repeating the analysis 
using hadrons in place of leptons and then scaling the yield by the momentum-dependent 
electron and muon misidentification probabilities. Table | summarizes the data yield and 
backgrounds. 

After subtracting the backgrounds, we are left with B — > D°Xlv and B — > D + X£v 
decays, where X stands for zero or more pions or photons. We divide each of these samples 
into ten equal bins of w in the range 1.0 < w < 1.6, where w is the reconstructed value 
of w, and is smeared by the detector resolution and motion of the B. In each w bin we 
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FIG. 1. The cos Q B -Di distribution for (a) B -> D°X£u and (b) B -> D+X£u. The data 
(solid circles) are overlaid with simulated B — > Dlu decays (dashed histogram), B — > decays 
(dotted histogram), 5 — > D**lv + D^iriu decays (dash-dotted histogram), and their total (solid 
histogram). The normalizations of the simulated samples are provided by the fit. 



fit the cos 6 b-di distributions of the B — ► D^Xiv and B — > D + Xlv samples for the yields 
of D£p, D*tv and the sum of D**iv and D^hiu, using Monte Carlo-simulated cos9 B _ D g 
distributions for each of these modes. In the simulation, we model B — ► .D£i/ decays using 
ISGW2 || and B — > £>*£z/ decays using the form factors measured by CLEO RlO| . We 
model the D** mesons with radial and angular excitations using ISGW2 and non-resonant 
D^'n states using the results of Goity and Roberts []Tlj| . The detector simulation is based 



on GEANT [|T2]]. In our fits, we apply the isospin symmetry constraint that the ratio of 
B — > D*iv to S — > -D£^ decay rates should be the same for charged and neutral S-meson 
decays. Figure [I] shows the result of a fit over all w. The sum of the B~ — > D°£~i? and 
B° — > D + £~v yields as a function of u> are displayed in Figure [| These w distributions are 
the basis of our studies of the B~ — > D°£~i? and 5° — ► D + £~i? form factor and decay rate. 
The differential decay width of B — > Div is given by M 



dV_ _ Gjr\Vc 12 

dw ~ 487T 3 



Fl " cb| (m B + mnfmUw 2 - l) 3/2 F D (wf, (2) 



where Gf is the weak coupling constant, mn is the mass of the D° or D + and Fd(w) is 
the form factor. We fit the B~ — > D°£~u and 5° — ► D + £rv yields in the intervals with 
to > 1.12 to extract information on |V^,| and the form factor. The two bins with w < 1.12 
are excluded because they suffer from small rates due to the (w 2 — 1) 3//2 suppression and 
large backgrounds. In our fit, the \ 2 function for B~ — > D°£~~u is expressed as 

10 \ATobs _ y-10 .../V.l 2 

where jV° 6s is the yield in the ith w) bin and 

N 3 = 2f + _N T(4S) B Kn T B - [ dwdT/dw (4) 
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FIG. 2. The sum of B~ — ► D®l~v and B° — > D + i~D yields as a function of io, for the data 
(solid circles) and using the best fit linear form factor (dashed histogram) or dispersion relation 
inspired form factor of Boyd et al. (solid histogram). 



is the number of decays in the jth w bin implied by the fit parameters. Here tb- is the B~ 
lifetime [13|, Bk-k is the D° — > K~n + branching fraction |l4j], A r x(45) is the number of T(45) 
events in the sample, and /+_ is the T(4S) — > B + B~ branching fraction. An efficiency 
matrix, ey, accounts for reconstruction efficiency and for the smearing of w. The fraction of 
decays in each w bin that are reconstructed ranges between 17% and 21% and the average 
w resolution is 0.026, about one-half the bin width. The small Monte Carlo statistical 
uncertainty in the efficiency matrix is represented by . We form Xo+e-p analogously. 
In the fits, we minimize \ 2 — Xd+£~p + xjjoi-py varying |T4&|-F£>(1), the coefficients in the 
parametrization of Fd(w) / Fo(l) , and /+_. We assume that the form factor parameters are 
common to D°£~i? and D + £~u decays and that B + B~ and B°B° together saturate T(4S) 
decays. 

We investigate several parametrizations of the form factor. Results for all of these are 
summarized in Table In all fits we find /+_ = 0.49 ± 0.04, consistent with previous 
measurements ||15|| , and similar correlation coefficients. We first consider the common ex- 
pansion F D {w) I F D {1) = 1 — p 2 D {w — 1) + Cb(w — l) 2 . When cd is constrained to be zero, 
we find p 2 D = 0.76 ± 0.16 and \V cb \F D {\) = 0.0405 ± 0.0045 with the correlation coefficients 
p(\V cb \F D (l),p 2 D ) = 0.95, p(\V cb \F D (l)J + ^) = 0.12 and p(f+-,p 2 D ) = 0.03. The X 2 is 8.8 
for 13 degrees of freedom. When the D°£~u and D + l~v samples are fit separately, they 
give consistent results for all parameters. This form factor is superimposed on the data in 
Figure When is allowed to vary, we find that it is consistent with zero within large 
errors; that is, our data allow substantial curvature but do not require it. In this fit, p 2 D and 
cd are completely correlated because our data are most precise at large values of w. 

Dispersion relations constrain the form factor. Boyd et al. |16j expand the form factor 
in the variable z = (y/w + 1 — \ / 2N)/(\/w + 1 + V2N), where iV ta 1.1. Because z is small, 
this expansion converges more rapidly than one in w — 1. Fitting for the linear coefficient a% 
with iV = 1.108, we find a\ = —0.043 ± 0.027. Expanding this form factor in powers of w — 1 
yields p 2 D = 1.30±0.27 and cd = 1.21±0.31 plus higher order terms. Caprini et al. [|ljj have 
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TABLE II. Summary of the B — * Dlv form factor fits. In addition to the quoted statisti- 
cal uncertainties, there are fractional systematic uncertainties of 12% in p D and cd and 8% in 
\V cb \F D {\). 



Form factor 


Pi 


CD 


W 2 \V cb \F D (l) 


X 2 /dof 


Linear 


0.76 ±0.16 




4.05 ± 0.45 


8.8/13 


Parabolic 


n 77+1-18 
u - ' ' -2.83 


01 + 1 - 70 


4 05+ 1 ' 51 


8.8/12 


Boyd et al.* ] 


1.30 ±0.27 


1.21 ±0.31 


4.48 ±0.61 


8.9/13 


Caprini et aU 


1.27 ±0.25 


1.18 ±0.26 


4.44 ±0.58 


8.9/13 



* We find a x = -0.043 ± 0.027 for N = 1.108. 



' This form factor also has terms of order (w — l) 3 and higher. 



also used dispersion relations to constrain the form factors. Their parametrization leads to 
similar results. 

To obtain the B —>■ Div decay rate, we use the form factor parameters provided by the fits 
and integrate dT/dw over to. The form factor of Boyd et al. [|16|] gives V = (14.1 ± 1.0) ns _1 . 
The other parametrizations give the same result within 1%. 

The systematic uncertainties are given in Table |T|. The uncertainties in the 5-meson 
momentum and mass dominate because they respectively affect the width and mean of 
the cos 9b-d£ distributions and therefore the D°fru and D + £~v yields extracted in each w 
bin. We have tuned our simulation to reproduce the B momentum distribution observed in 
fully reconstructed B decays; however, there is a 6 MeV uncertainty in the mean and this 
leads to fractional systematic uncertainties of 5% for p 2 D , 4% for |T^,|Fd(1) and 3% for T. 
The 1.8 MeV [HJ uncertainty in the B mass generates uncertainties of 7% for p 2 D , 4% for 
1^1^(1) and 3% for V. 

The other large systematic error arises from uncertainty in the cos 9 b-di distribution of 
the combined B — > D**£u and B — > D^7r£u backgrounds. This distribution depends mainly 
on the number of final-state pions that are not reconstructed. We therefore separate it into 
two components: one in which the final-state D is accompanied by one 7r and the other in 
which it is accompanied by two. We vary their relative proportions from 1:4 to 2:3 [19] to 
evaluate the systematic uncertainty. 

The B — > D*£u form factors affect the distribution of these decays in cos6b~di and 
therefore influence the extracted B — > Div yield. We vary the form factors within the 
uncertainties of the CLEO measurement flO|| , taking into account the correlations among 
the form factor parameters (i?i,i?2 and p 2 ). 

The slope of the B —>■ Dlv form factor is varied in our fit so its uncertainty is included 
in the statistical error in the decay width. Using a linear form factor to extract the width 
rather than the dispersion-relation-inspired form factor changes the width by 0.8%. The 
form factor can also affect the cos 6 ] b-di distribution used to extract the D®l~v and D + £~i? 
yields in each w bin and the efficiency matrix. Each of these effects is less than 1%. 

Our final result is 

F(B -> Dlv) = (14.1 ± 1.0 ± 1.2) mT 1 . (5) 
Multiplying this by the measured 5-meson lifetimes gives the branching fractions 
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TABLE III. The fractional systematic uncertainties. 



Source 




l^l-Pb(i) 


Y{B -> Dlv) 


Track-finding 




0.02 


0.035 


Lepton ID 


— 


0.01 


0.020 


K and tt ID 


0.02 


0.01 


0.022 


Backgrounds 


0.06 


0.04 


0.018 


\pb\ and Mb 


0.08 


0.05 


0.042 


Luminosity 


_ 


0.01 


0.018 


Dlv form factor 


0.01 


0.01 


0.010 


D*lv form factors 


0.01 


0.01 


0.005 


D**£u model 


0.04 


0.03 


0.026 


D branching fractions 




0.02 


0.036 


TB 




0.02 


0.026 


Total 


0.11 


0.08 


0.085 



B(B~ -> D°£-9) = (2.32 ± 0.17 ± 0.20)% and (6) 
B(B° -> D + t~v) = (2.20 ± 0.16 ± 0.19)%, (7) 

where the first errors are statistical and the second are systematic. Since we derive both 
branching fractions from the decay width, their errors are completely correlated. This result 
is consistent with previous measurements but is more precise. Combining it with the previous 
CLEO measurement B, taking into account statistical and systematic correlations, gives 

T(B -> Dlv) = (13.4 ±0.8 ± 1.2) ns _1 , (8) 
B(B~ -»• D°£~v) = (2.21 ±0.13 ±0.19)%, and (9) 
B(B° -> D + l~v) = (2.09 ± 0.13 ± 0.18)%, (10) 

where the errors in the partial width and branching fractions are completely correlated. 

Our studies of the form factor give p 2 D = 0.76 ± 0.16 ± 0.08 (linear fit) and p 2 D = 1.30 ± 
0.27 ± 0.14 and cd = 1.21 ± 0.31 ± 0.15 plus higher order terms (dispersion relations). The 
latter gives 

\Vcb\F D (l) = (4.48 ± 0.61 ± 0.37) x 10~ 2 . (11) 

Various authors have found F D (1) = 0.98±0.07 and F D (1) = 1.03±0.07 [|, and a recent 
lattice calculation finds the preliminary value -Fd(I) = 1.069±0.029 0]. Using ^(1) = 1.0, 
we find \V c b\ = 0.045 ± 0.006 ± 0.004 ± 0.005, where the last uncertainty covers all of these 
values of Fd(1). This value of \V c b\ is consistent with that from B — > D*lv decays, though 
its uncertainty is larger. Using the linear form factor, as has been done in most previous 
studies of B — > D*£v, gives a value of \V c b\ that is about 10% smaller than this. While the 
curvature of the form factor is likely to have a smaller effect on the \V c b\ extracted from 
B — > D*iv decays, its effect could nevertheless be important. 
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